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Twofold spontaneous symmetry breaking in a heavy fermion superconductor UPt 3 
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The field-orientation dependent thermal conductivity of the heavy- fermion superconductor UPt3 
was measured down to very low temperatures and under magnetic fields throughout three distinct 
superconducting phases: A, B, and C phases. In the C phase, a striking twofold oscillation of the 
thermal conductivity within the basal plane is resolved reflecting the superconducting gap structure 
with a line of node along the a axis. Moreover, we find an abrupt vanishing of the oscillation across 
a transition to the B phase, as a clear indication of a change of gap symmetries. We also identify 
extra two line nodes below and above the equator in both B and C phases. From these results 
together with the symmetry consideration, the gap function of UPt3 is conclusively determined as a 
E\ u representation characterized by a combination of two line nodes at the tropics and point nodes 
at the poles. 
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Spontaneous symmetry breaking is one of the funda- 
mental paradigms encompassing from condensed matter 
physics to high energy physics, constituting the founda- 
tion of modern physics. This paradigm is crucial some- 
times because it can give a handle to discover some un- 
known exotic ordered phase. This is particularly true 
when broken symmetry is extremely low, that is, the 
"residual symmetry" is so small, one may effectively and 
self-evidently narrows down possible ordered phase to 
identify. 

Understanding the unconventional superconductivity, 
in which electron pairs are formed without phonon, has 
been a challenge. Part of the problem in uncovering the 
mechanism is that little is known about the pairing sym- 
metry. The heavy-fermion superconductor UPt3 is one 
of the examples whose pairing symmetries are as yet to 
be clarified. The most intriguing feature of this material 
is the existence of a multiple phase diagram; UPt3 un- 
dergoes a double superconducting transition at the upper 
critical temperature T c + ~ 540 mK into the A phase and 
at the lower critical temperature T~ ~ 490 mK into the 
B phase [1]. In addition, the third (C) phase is stabi- 
lized at low temperatures under high magnetic fields [2]. 
A crucial role of a weak antiferromagnetic order below 
Tn ~ 5 K for thephase multiplicity is indicated by the 
pressure studies [3[ . Power law dependence of the ther- 
modynamic and transport quantities reveal the presence 
of nodes in the superconducting gap [4H7J. Moreover, a 
possibility of an odd-parity pairing is inferred from the 
nuclear magnetic resonance studies of the Knight shift [8] 
and is supported theoretically [9] by eliminating the sin- 
glet even parity scenario. 

Extensive theoretical efforts have been devoted to ex- 
plain these disparate experimental results [9l4l l| . Among 



them, the E^ u scenario with a line node in the basal 
plane and point nodes along the c axis has been regarded 
as one of the promising candidates [12j. Several exper- 
imental results, such as the anisotropy of the thermal 
conductivity [13| and the ultrasonic attenuatio n | 14J| as 
well as the recent small- angle neutron scattering [15| and 
the Josephson tunnel junction [16j, have been claimed 
to be compatible with this model. On the other hand, 
there exist some controversies in explaining the following 
experiments; 1) the spontaneous internal field due to the 
broken time-reversal symmetry is most likely absen t) 17j , 
2) the d- vector has two components in the B phase [8|, 3) 
a point where the three superconducting phases meet is a 
tetracritical point [2]. Moreover, to date no experimental 
evidence for the gap structure of each phase associated 
with the E2u model has been provided. The pairing sym- 
metry of UPt3, therefore, remains unclear. 

One of the most conclusive ways to identify the pairing 
symmetry is to elucidate the gap structure by the ther- 
mal conductivity measurements with rotating magnetic 
fields relative to the crystal axes deep inside the super- 
conducting state. This technique has been successful to 
probe the nodal gap structure of several unconventional 
superconductors by virtue of its directional nature and 
sensitivity to the delocalized quasiparticles [18|. In this 
paper, we present a decisive experiment of the angular 
dependence of the thermal conductivity of UPt3 reveal- 
ing the spontaneous rotation symmetry lowering, namely 
the unusual gap structure with a lower rotational sym- 
metry than the crystal structure. 

High quality single crystal of UPt3 with the high resid- 
ual resistivity ratio of 800 was grown by the Czochralski 
pulling method in a tetra-arc furnace [19|. We measured 
the thermal conductivity along the hexagonal c axis (heat 



current q \\ c) on the sample with a rectangular shape (3 
x 0.42 x 0.4 mm 3 ). To apply the magnetic fields with 
high accuracy relative to the crystal axes, we used a sys- 
tem with two superconducting magnets generating the 
fields in two mutually orthogonal directions. The mag- 
nets are installed in a Dewar seating on a mechanical 
rotating stage, enabling the continuous rotation of the 
magnetic fields. 

First, we begin with demonstrating that the thermal 
conductivity (n) well probes the superconducting quasi- 
particle (QP) structures from its temperature (T) and 
magnetic field (H) dependences. From now on, the 
hexagonal [1210], [1010], and [0001] axes are denoted as 
the a, 6, and c axes, respectively. The inset of Fig. 1 
shows the T dependence of n{T)/T under zero field and 
3 T along the b axis. With decreasing T, the zero-field 
k(T)/T shows a steep increase up to ~ 0.3 K without 
apparent anomalies at T c + and T~ . On further cool- 
ing, k(T)/T considerably decreases due to a reduction 
of the QP densities, and takes an extremely small value 
at the lowest T ~ T+/20, consistent with the previous 
measurements [7J]. In the normal state (3T), k(T)/T ap- 
pears to continuously increase down to the lowest T. The 
dashed line denotes n(T)/T obtained from the normal- 
state resistivity p(T) using the Wiedemann-Franz law, 
k(T)/T = Lq/p(T) (Lq\ the Lorentz number). Impor- 
tantly, we confirm that n(T)/T is close to Lo/p(T) at 
low temperature T < 100 mK, indicating the dominant 
electronic contribution in the heat transport. In this 
T-range, the H dependence of the thermal conductivity 
k{H)/T at 55 mK shows a remarkable i^-linear depen- 
dence at low fields for both c and b directions (the main 
panel of Fig. 1) in contradiction to the field-insensitive 
behavior of fully gapped superconductors except in the 
vicinity of H C 2 [20[ , providing evidence for the nodal su- 
perconductivity in UPt3. 

In addition, we find distinct anomalies associated with 
a transition from the B to C phase at Hbc (open arrows). 
The fact that the BC transition manifests by a sharp 
change of the slope implies a suppression of one of the 
degenerate order parameter components in the B phase. 
This behavior can be more clearly resolved for the b axis. 
The determined Hbc together with H C 2 denoted by the 
solid arrows are summarized in Fig. 3(d) for H \\ b. We 
also note that a striking anisotropy is found in k(H)/T 
at 55 mK near H C 2'. n/T forH || c shows a rapid increase 
just below iJ C 2, while the one for H || b linearly increases 
up to H C 2, as similarly observed in S^RuC^ [21(. A 
search of the relevance of this behavior to the odd-parity 
superconductivity is a fascinating issue to be addressed. 

Next, to shed light on the nodal topology in the su- 
perconducting phases, we concentrate on the angular de- 
pendence of k. The most significant effect on the thermal 
transport for nodal superconductors in the mixed state 
comes from the Doppler shift of the QP energy spectrum, 
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FIG. 1: (color online). Magnetic field dependence of the ther- 
mal conductivity k(H)/T along the c and b axes at various 
temperatures. The open and closed arrows represent the B — >• 
C transitions Hbc and the upper critical fields if C 2, respec- 
tively. Inset: temperature dependence of k,(T)/T under zero 
field and at 3 T for H || b. The dashed line shows k/T = Lo/p 
(Lo: the Lorentz number) obtained from the normal-state re- 
sistivity p using the Wiedemann-Franz law. 



E(p) —> E(p) — v s • p, in the circulating supercurrent 
flow v 3 . This effect becomes important at such posi- 
tions where the gap becomes smaller than the Doppler 
shift term (A < v a • p). The maximal magnitude of the 
Doppler shift strongly depends on the angle between the 
node direction and i7, giving rise to the oscillation of 
the density of states (DOS). Consequently, n attains the 
maximum (minimum) value when H is directed to the 
antinodal (nodal) directions [22j. Figure 2 shows n{(j)) 
normalized by the normal state value n n as a function 
of the azimuthal angle <p at 50 mK (~ T+/10) at \p^H\ 
= (a) 3.0 T, (b) 1.0 T, and (c) 0.5 T, respectively. The 
data are taken in rotating H after field cooling at <p = 
-70°, and n n is measured at 50 mK above H C 2 for H \\ b. 
In the normal state (3.0 T) and the B phase (0.5 T), we 
find no ^-dependence within experimental error. 

By contrast, what is remarkably is that k((/)) exhibits a 
distinct twofold oscillation with a minimum at <\> = 0° in 
the C phase (1.0 T). The open circles are obtained under 
field cooling condition at each angle. The data obtained 
by different procedures of field cooling coincide well with 
each other, indicating negligibly small effect of the vor- 
tex pinning. Strikingly, since the twofold symmetry is 
lower than the hexagonal crystal structure^ the in-plane 
anisotropy of the Fermi surface and H C 2 [l2[ is immedi- 
ately ruled out as the origin. As shown by the solid lines, 
K,((f)) can be decomposed into two terms; k((/)) = fto + ^20, 
where ^o is a ^-independent term and ^20 = C^ cos 20 
is a twofold component. Figure 2(e) shows the ampli- 
tude of the twofold component \C2<f>/^n\ as a function of 
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FIG. 2: (color online). Angular variation of the thermal con- 
ductivity ft(</>) normalized by tt n at 50 mK as a function of the 
azimuthal angle <j> for \/jloH\ = (a) 3.0 T, (b) 1.0 T, and (c) 
0.5 T, respectively. «(0) is measured with rotating H within 
the ab plane (the polar angle = 90°) as schematically shown 
in (d), where cj) and are measured from the a and c axes, 
respectively, and q is injected along the c axis. The solid lines 
show the twofold component in «(0)/« n . The open circles 
represent «(0)/« n at 1 T obtained under the field cooling 
condition at every angle, (e) Field variation of the twofold 
amplitude |C2</>/ftn| at 50 mK at — 90° (solid circles) and 
63° (open circle), respectively. 



H/H C 2, where H C 2 = 2.6 T for H || b. It can be clearly 
seen that IC^/ftnl suddenly appears to be finite ~ 3% 
in the C phase, implying a change of the gap symmetries 
across the BC transition that is of second order. We note 
that \C2cf)/^n\ obtained by rotating H conically around 
the c axis at fixed = 63° is same order of magnitude 
with the values at = 90° as denoted by an open circle 
in Fig. 2(e). 

To further elucidate the gap symmetry, we present 
the polar angle (0) dependence of n in Fig. 3, showing 
K,(6)/K n measured by rotating H within the ac plane 
(green circles) and the be plane (orange circles) at 50 
mK at \fji H\ = (a) 1.5 T, (b) 1.0 T, and (c) 0.5 T. Here, 
n n is measured at 50 mK above H C 2 for H || c. The dom- 
inant twofold oscillation is found in all the fields with 
maxima at = 90°, which could be attribute to, such as 
the Fermi surface and/or the gap topology or the differ- 
ence in transport with H parallel to and normal to the 
heat current q. Regardless of the origin, the fact that 
n(0)/ n n is maximized at = 90° excludes an artificial 
origin of the in-plane twofold oscillation in the C phase 
due to a misalignment of H relative to q. We thus con- 
clude that the in-plane twofold symmetry in the C phase 
is a consequence of the node. 

In the B phase (0.5 T), the two different scanning pro- 
cedures within the ac and be planes well converge with 
each other, consistent with the ^-independence of k. In 
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FIG. 3: (color online). Angular variation of the thermal con- 
ductivity k(6) normalized by K n at 50 mK as a function of 
the polar angle for \fi H\ = (a) 1.5 T, (b) 1.0 T, and (c) 
0.5 T, respectively. The K,(0)/K n curves measured by rotating 
H (inset of (a)) within the ac plane (green circles) and the 
be plane (orange circles) are simultaneously plotted. Inset of 
(c): Aft(#)/ft n = (k(0) — ko — K2o)/Kn vs plot at 50 mK at 
0.5 T, where ko is a ^-independent term and K20 = C26 cos 20 
is a twofold component . (d) The phase diagram of UPt3 
with the three superconducting phases, labelled A, B, and C, 
for H || b. The red and blue circles represent Hbc and H C 2, 
respectively, deduced from the present measurements. The 
schematic shapes of the gap symmetries for each phase are 
shown. 



addition, we find extra two minima at = 20° and 160°. 
By plotting Aft(#)/ft n = (k(0) — fto — ^26>)/^n vs 6 after 
the subtraction of ^0 and K2Q = C2Q cos 20, the minima 
become clearly visible at 35° and 155° (Fig. 3(c), inset). 
This double-minimum structure is also found in the C 
phase (Fig. 3(a)). We infer that these minima are de- 
rived from the two horizontal line nodes at the tropics as 
discussed below. In contrast to the B phase, the two scan- 
ning results do not coincide in the C phase (Fig. 3(a)); the 
difference is diminished at the poles and maximized at 
= 90°, being consistent with the in-plane twofold sym- 
metry. Moreover, a significant appearance of the twofold 
symmetry across the BC transition can be seen at 1.0 T 
(Fig. 3(b)), in which one experiences the BC (CB) transi- 
tion twice by varying because of the anisotropy of Hbc • 
Indeed, the transitions occur at = 30° and 150° taking 
distinct kinks. Remarkably, the difference between the 
two scanning procedures becomes finite upon entering 
the C phase, providing the compelling evidence for the 
twofold symmetry of the gap structure in the C phase. 
Moreover, the fact that |C2</>/ft n | takes same order of the 
magnitude at = 90° and 63° is in favor of a line node 
along the a axis rather than the point nodes in the basal 
plane. Notably, although a mechanism which fixes do- 
mains is a puzzle, the in-plane twofold symmetry of n{(j)) 
indicates a single superconducting domain. 



We discuss the order parameter symmetry of UPt3 
within the triplet category. The present experiments in- 
dicate (i) the line node along the a axis in the C phase, 
(ii) the absence of in-plane gap anisotropy in the B phase, 
and (hi) the two line nodes at the tropics in both B and 
C phases. Taking into account all these results and the 
d- vector configurations assigned by the Knight shift [8|, 
the order parameter is unambiguously determined with 
a form of (k a b + kbc)(5kl — 1) for the B phase, where b 
and c are unit vectors of the hexagonal axes representing 
the directions of d- vectors. This state belongs to two- 
dimensional E\ u representation with the /-wave charac- 
ter, the so-called planar state in triplet pairing in the D§h 
hexagonal symmetry, and to degenerate E u state for the 



recent claimed D^d trigonal symmetry [23|, |24[ . The gap 
structure consists of the two horizontal line nodes at the 
tropics (k c = ±l/y/E, = 63° and 117°) and the point 
nodes at the poles (k a = fc& = 0). Note that although 
the locations of the horizontal line nodes estimated by 
assuming a spherical Fermi surface do not agree with the 
observation (0 = 35° and 155°), it could be changed by 
considering the realistic Fermi surface [12J . 

By lifting the doubly degeneracy, the order parameter 
for the C phase is given by kbd(5k^ — 1) for H \\ ab and 
k}jd(bk^ — 1) for H || c, respectively. In the same manner, 
k a b(JSk1 — 1) state is readily assigned for the A phase. 
The schematic shapes of the gap symmetries in the three 
phases are shown in Fig. 3(d). We emphasize that this 
state is compatible not only with the hybrid gap state 
indicated by the several experiments [13|, LL4| , in the sense 
that the line and point nodes simultaneously exist, but 
also with some experimental results for which the E2 U 
model [9] has failed to describe, i.e., the absence of the 
internal magnetic field [17j, the two-component d- vector 
for the B phase [8|, and the tetracritical point in the 
phase diagram [2[. 

To further strengthen our identification, in particular 
on the existence of the horizontal line nodes on the trop- 
ics, we calculate the angle-resolved DOS by solving the 
Eilenberger equation [25] for several possible gap func- 
tions. We compare here putative three gap functions in 
the C phase relative to the data in Fig. 4 where k(0)/k u 
and the DOS differences along the vertical nodal and 
antinodal ^-scannings are depicted. The double peak 
structure characteristic in E2 U and E\ g whose origin 
comes from the horizontal node on the equator is not sup- 
ported by the data that are consistent with the present 
E\ u with the horizontal nodes on the tropics. In view of 
the Doppler shift idea mentioned above the QPs in the 
horizontal node on the equator contribute more when the 
field direction is away from = 90°. 

In summary, we find striking twofold oscillations in 
angle-resolved thermal conductivity measurements at low 
temperatures in a strongly correlated heavy fermion su- 
perconductor UPt3. This spontaneous symmetry low- 
ering, which is the lowest possible rotational symmetry 
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FIG. 4: (color online), ^-dependence of the thermal conduc- 
tivity obtained by subtracting the green data from the orange 
data in Fig. 3(a) (left axis) and the density of states difference 
normalized at = 90° (right axis, arbitrary scale) along the 
vertical nodal and antinodal scannings for three possible gap 
functions in the C phase: 1. The present E\ u (kb(5kc — 1)), 
2. E\ g (kbk c ), 3. E2u (kakbkc). Those gap structures are 
sketched in the inset. 



breaking in hexagonal crystal fortuitously and effectively 
narrows down the possible symmetry classes and leads us 
to uniquely identify the pairing symmetry for each phase 
in the multiple phase diagram. We conclude that the 
realized pairing function is E\ u with the /-wave charac- 
ter, i.e., the so-called planar state in the triplet pairing. 
This state is analogous to the B phase in superfluid 3 He, 
and obviously bears the Majorana zero mode at a sur- 
face [2a, l22| , namely a topological superconductor that 
is quite rare to find. Thus it is worth exploring further 
to understand this interesting material as a new platform 
for topological physics. 
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